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Abstract. We shall define the g-analogs of multiple zeta functions and multiple polylogarithms 
in this paper and study their properties, based on the work of Kaneko et al. and Schlesinger, 
respectively. 

1 Introduction and definitions 

Let < g < 1 and for any positive integer k define its g-analog [k] — [k]q — (1 — q'')/{l — q). In 
[S] Kaneko et al. define a function of two complex variables fq{s; t) = J^kLi '?'°V[^]* such that the 
^-analog of Riemann zeta function is realized as 

Cq{s) ■■= fq{s;S- 1). 

For any property V let be the set of integers satisfying V. After analytically continuing fq{s; t) 
to as a meromorphic function Kaneko et al. proved the following main result 

Theorem 1.1. (0 Prop. 2, Thm. 2]) One can analytically continue Cq{s) to C \ {Z<o + i^|^^#o}- 
Moreover, for any s 7^ 1 one has 

limCqis) ^ a^). 

They also study the special values of Cq{s) at non- negative integers. In this paper we shall gener- 
alize these to the (Euler-Zagier) multiple zeta functions, which are defined as nested generalizations 
of Riemann zeta function 

C(si,...,s,)= k^'-'-kT' (1) 

0<fei<---<fcd 

for complex variables Si, . . . , s,; satisfying + • • • + ct^ > d — j + \ for all j — 1, . . . , c?. Here 
and in what follows, whenever s G C we always write a — 5R(s), the real part of s. The analytic 
continuation of multiple zeta functions has been studied independently in (Jj and We know 
that C(sij ■ • • J Sd) can be extended to a meromorphic function on \ 6^ where 



6<j = < (si, . . . , Sd) G 



Sd =1, Sd-i + Sd = 2, 1, -2m Vto e Z>o, 
and + • • • + G Z<d_j+i < d — 2. 



(2) 



To find the g-analog of multiple zeta functions we first define an auxiliary function of 2d complex 
variables si, . . . , Sd,ti, . . . ,td G C 

Eqkiti^ hkdtd 
r, — imr 

which converges if ^{tj + • • • + i^) > for all j = 1, . . . , (see Prop. 12^ . In the next section we 
are going to analytically continue this function to C^"^ as a meromorphic function with explicitly 
defined poles. 
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We now define the g- multiple zeta function by specialization of fq: 

Cq(si, . . . , Sd) := fqi-Sl, . . . , Sd;Sl - 1, . . . ,Sd - 1) 

which will be shown to be the correct g-analog of multiple zeta functions. When aj > 1 for all j we 
can express this by the series 

qki(si-l)-\ hfcd(sd-l) 

0</ci<---<fcd ' ' 

Note that when d — 1 this is the same as the g-analog of the Riemann zeta function defined in 
Put 

Sd e 1 + Z, or Sd e Z<o + Z^o, 

ft' J f \ ^(Td log*? " log*? I ^ « 

or H h Sd e Z<d-j+i + Z, j < 

logg j 

Here the last part in 6^ is vacuous li d ~ 1. The primary goal of this paper is to prove 

Main Theorem. The q-multiple zeta function Q(si, . . . , Sd) can be extended to a meromorphic 
function on \ S'^ with simple poles along S^. Further, for all {si, . . . , Sd) £ C"^ \ &d 

limCg(si,...,Sd) = C(si,...,Sd)- 

9T1 



In sectional we propose a new definition of the g-multiple polylogarithms and briefly study their 
properties. We also review Jackson's g-derivatives and g-definite integrals and define g-iterated 
integrals as g-analogs of Chen's iterated integrals. 

It is known that there're two kinds of shuffle relations among multiple zeta values (MZV for 
short). The first one is produced by their power series expansions, the second by using Chen's 
iterated integrals. In the last section of this paper we will apply our g-iterated integral technique 
to g-multiple polylogarithms in order to study the g-shuffle relations of the second kind for g-MZV. 
For simplicity we will only deal with C,q[m)C,q{n) for positive integers m ^ n. These relations reduce 
to the ordinary ones when g t 1- We thank Prof. Kaneko for his questions relating to this part of 
our study and sending us his offprint upon which the current work is based. 



2 Analytic continuations of fq and Q 

The purpose of this section is two-fold: first we will use the auxiliary functions fq introduced in the 
first section to give a quick analytic continuation of g-multiple zeta functions Cg(si, . . . , Sd), though 
this is not enough to show it's the right g-analog of the multiple zeta functions. Second, we write 
down these expressions involving binomial coefficients explicitly which will be used to study special 
values of C(j(si: ■ • ■ j Sd) in section|Hl 
We need a simple lemma first. 

Lemma 2.1. Let k he a positive integer. For all 1 > q > 1/2 we have 

1 - g'= 

1 < ^ < 2k. 

1 - g 

Proof The inequality 1 < is obvious. Let /(g) = 2fc(l - g) - (1 - g'^). We show that /(g) > 
for all < g < 1. Indeed, for all such g 

/'(g) = -2k + kq^-^ < -fc < 0. 

So / is strictly decreasing and the positivity of / follows from /(I) = 0. This implies that 2k > 
(1 - g*')/(l - g) as desired. □ 
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Proposition 2.2. The function 



^fcltlH hfcrftd 



0<fel<---</£d 



converges if^{tj + - ■ ■+tj) > for all j = 1, . . . , d. It can he analytically continued to a meromorphic 
function over C'^'^ via the series expansion 



fq{si,...,Sd;ti,...,td) = (l-q) 



wt(s ) 



+CXD d r 

E n 

ri,...,rd=0 j = l 



Sj + - 1 



1 — q^:i+*i^ \-rd+td 



(4) 



where wt(s ) — si + ■ ■ ■ + Sd- It has the following (simple) poles: tj + --- + tdG Z<o + y^^^ 
Proof. Assume |5R(sj)| < Nj and let t, = ^{tj) for all j = 1, ... , d. By Lemma 



f,isi,...,sd;ti,...,td) < J2 llC^k.f'l''"' 

0<ki<-<kdj=l 



(5) 



Let k — kd-i (when d — 1 take kg — 0), n — kd, N = Nd, and t = Td- Then by root test 
J2n>k iT'^Q^^ converges and moreover 

N 

" "• ■ (6) 



n>k ^ ^ n>k 



where f]y{x;q,T) = 'l2iLo'^i^^ ^ polynomial of degree N whose coefficients depend only on the 
constants N, q and r. This proves the first part of the lemma when d = 1. In the general case it 
follows from (|SJ) and © that 



fq{si,...,Sd;ti,...,td) < Q _ r\N E'^' 



1=0 0<ki<---<kd-i<k \j = l 



Hence the first part of the lemma follows from an easy induction on d. 
By binomial expansion (1 — x)^^ — X]^o i^^r^^)^^ S"^^ 



+ 00 



f,is^,...,Sd;t,,...,td) = {l-qr<^'^ E H 

0<ki<---<kdru---,rd=0 \j = l 



S,- + r,- - 1 



ri 



As < 5 < 1 the series converges absolutely by Stirling's formula so we can exchange the sum- 
mations. The proposition follows immediately from the next lemma by taking Xj = for 

Lemma 2.3. Let Xj G C such that \xj \ < 1 for j — 1, . . . , d. Then 



0<ki< - <kd j=l 



3 = i 



'^j ' ' ' 

J ■ ■ ■ Xd 



(7) 



Proof. When n = 1 this is clear. For n > 2 we denote the right hand side of {Tj) by Fd{xi, . . . , Xd). 
Then 



d-l 



E m= E E 

0<ki<---<kd 3 = 1 0<ki<---<kd-i 3 = 1 kd>kd-i 



Xd 



Xd 



-Fd-lixi, . ■ . ,Xd-2,Xd-lXd)- 



The lemma follows immediately by induction. 



□ 
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Recall that CgC'Si, ■ ■ ■ i Sd) ~ fqi^i, ■ • ■ , Sd; si — 1, . . . , Sd — 1). Hence we have the following 
immediate consequence. 

Theorem 2.4. The q-multiple zeta function Cg(si, . . . , Sd) can be extended to a meromorphic func- 
tion with simple poles lying along 6^; 



c,{su...,sd) = {i^qr^'^ J2 n 

ri....,rd=Q j=l 



1 — q^3+^j^ hi-d+Sd-d+j-l 



To see the effect of taking different specializations of tj in /g we define the shifting operators Sj 
(1 < j < d) on the multiple zeta functions as follows: 

<SjC{si, ...,Sd)= C(si, . . . ,Sd) + {I ~ q)C{si, . . . , - 1, . . . , Sd). (8) 

It is obvious that these operators are commutative. 

Proposition 2.5. Let ni, . . . ,nd be non-negative integers. Then we have 

fq{si, . . . , Sd] Sl - 1 ~ Hi, . . . , Sd - 1 ~ Ud) 

Jii rid / d 

r 



S^^o...oS2%isu...,Sd)^J2---J2 \ I[(j.)(^-iy' I Qisi- n,...,sd~rd). 

ri=0 rd=0 \J = 1 



Proof. We only sketch the proof in the case ni = ■ ■ ■ = Ud-i = 0. The general case is completely 
similar. In the rest of the paper we always let S be the shifting operator on the last variable. 
Suppose Ud = n = I. Then 

/g(si, . . . ,Sd; Si - 1, . . . , Sd-i - l,sd- 2) 

qki(si~l}-\ hfcd-i(sd_i-l)+fcd(sd-2) 



E 



qki{si — l)-\ h*;d-i(sd-i — 1) gfed(sd-2)^2 — q^'') + qkd{sd-i) 



=SC{si,...,Sd). 

The rest follows easily by induction. □ 

The next corollary answers an implicit question in jSj. 
Corollary 2.6. Let n be a positive integer. The specialization oft in fq{s;t) to s — I — n is 



fqis; s-l-n) = S-Cqis) = E ( J " 



We observe that one effect of the shifting operator is to bring in more poles. Essentially, 5" 
shifts all the poles of Cq{s) by n to the right on the complex plane. 

3 Analytic continuation of multiple zeta functions 

Let's begin with a review of some classical results on Bernoulli polynomials Bk(x) which are defined 
by the generating function 



fe=0 

Let Bk{x) be the "periodic Bernoulli polynomial" 

Bk{x) ^ Bk{{x}), x>l, 
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where {x} is the fractional part of x. Then we have Ch. IX, Misc. Ex. 12]) 

nez\{o} ^ ' 

Recall that the Bernoulli numbers satisfy i?^ = if fc > 2 while = 1 and Bi = —1/2. 

Lemma 3.1. For every positive integer M > 2 and x > 1 we have 

AMI 



Bm{x) 



< 



(27r)A^ ■ 

Proof. It follows from the fact that ({M) < C(2) = 7rV6 < 2 for M > 2. □ 

We know that one can analytically continue the multiple zeta functions as independently pre- 
sented in QjC] and [J by different methods. Moreover, ^{si, . . . , Sd) has singularities on the hyper- 
planes in 6d defined by (0). However, neither approach is suitable for our purpose here. So we 
follow the idea in ^ to provide a third approach in the rest of this section. The same idea will also 
be used to deal with the g-multiple zeta functions. 

Let's recall the classical Euler-Maclaurin summation formula TOI, 7.21]. Let f{x) be any 
(complex- valued) C°° function on [l,oo) and let m and M be two positive integers. Then we 
have 

n=l "'l r=l ^ '' 

/ BM+i{x)f^''+^\x)dx. (10) 



To simplify our notation, in definition we replace s^, ^d-i, kd by s, k and n, respectively. 
Taking f{x) = l/a;** and m = k and oo in H10[l we have: 



— ' n" ^ — ' n" 

n>k n=l n=l 



r fix) dx - ifik) y-^^f(r)(k) ~ ^rST.. \ 

Jk ^ ' 2 ^ (r + W (M + 1)! 



^ ' ' BM+iix)f'^''+'\x)dx 



J L I V ^'■+1 _ (g)M+l Z""" Sm+i(x) 



dx. 



{s-l)k''-^ 2k-' ' ^^{r + iy.k''+'- {M + iy.Jk x^+A'-f+i 
Here we have used the fact that Sfe = if A; > 3 is odd. By definition ^ we have 
Theorem 3.2. For all {si, . . . ,Sd) G \ &d and M > 1 + \ad\ + Wd-i\ we have 

M+l ^ 

C(S1, ...,Sd)= ^ —^{Sd)r-1 ■ C{S1, Sd~l + Sd + r -1) 



r\ 

r=0 



{sd)M+i 1 f°° Bm+i{x) 



where we set (s)o — 1 and (s)-i — l/(s — 1). This provides an analytic continuation of C,{si, . . . ,Sd) 

to c^\ed. 

Proof. We only need to show that the series in converges. Lemma f3.ll implies (if d = 2 then 
take fco = 1) 



E 



1 p Bm{x) , 
' dx 



4M! ^ 



(2^)M(Af-|a,|) ^ ^M-\.a-A-W.\ 



kd-i=kd-2 

which converges absolutely whenever M>l-|-|(T£(| + |o'rf_i|. □ 
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4 Proof of Main Theorem 

Fix (si, . . . , Sd) G such that Cj + ■ • • + 0"^ < d — j + 1 for all j — 1, . . . ,d. When d — 1 this 
is Thm. 11.11 due to Kaneko et al. We now assume d > 2 and proceed by induction. The key is 
a recursive formula for Cg(si, . . . , Sd) similar to Hll|) for ■ • ■ , Sd)- To derive this formula we 
appeal to the Euler-Maclaurin summation formula 1)11) |l again. Hence we set 

F{x) 



(1-g- 
as in 0. Then 

F'{x) =(logg)g^(^-i)-''"^ + '^" 



x\2 



In definition |j2Jl we replace Sd, kd-i, and kd by s, fc and n, respectively. We now take M = 1, 
f{x) — F{x + fc — 1) and let m ^ co in and get 

n>k ' \ n=l / 

= {1 - qy ^j^ F{x) dx ~ \F{k) - ^F'{k) B,{x)f"ix) dx) 

(g-1) igfc(^-i) 1 logg (7'=("-i)(s + g''' - 1) 



(s-l)logg [fc]''-! 2 [kY I2q-l [/c]^+i 

- " - f + " - " - (12) 

because B2{x+k — l) — B2{x) by periodicity. By the same argument as in |3], setting the incomplete 
beta integrals 

bt{a,l3) = I u'^-^il-uf-^du, 



we can obtain from ^ the following expression for the last term in (|12|l : 

_ ^ (1-griogg ^ a,(5)V(5- 1 + ^^,-5 + ^) (13) 

nGZ\{0} ^ ' v=±\,<d 

where (5 — 27rj/logg, a_i(s) = s(s+l),ao(s) = — 3s,andai(s) = 1. Repeatedly applying integration 
by parts on these incomplete beta integrals we get for v — ±1,0 and positive integer M >2 

hAs-\ + 8n,-s^y)= V (-1)^-^ + .7 .. ,^ ^ 
^ ^ (s - 1 + &n)r (1 - 9'=)«+'-'' 

^^^_l_(sJ_j__^OM-l 

(s - 1 + (5n)j\/-i 



(-l)^^- \\ ^ V (s - 2 + Af + -s - M + 1 + i^). (14) 



Set s' = (si, . . . , Sd-2) if d > 3 and s' = if d = 2. Putting ©, 11211, ^ and lO together 
and applying Prop. 1^51 we get 

Cg(s ', Sd-1, Sd) = , .J^} Cg(s", Sd-1 + Sd - 1) - ^5Cq(s', + Sd) 

(s-l)logg 2 
+ ^^52C,(s',s,_i + s, + l) + i^5C,(s',s,_i+s,)- (^-+^-) (15) 
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where Ci, and are the contributions from the sum involving bqk{. . . , —s + v). Exphcitly they 
are computed as foUows. Write 



r-l 



T{q,s,n,r) = Y[{2nin+{s-l+j)\ogq) \ (16) 



Then 



M-l ^ / , X r+1 

T{q,s,n,r) f logq \ , cS/- , , l1^ 

r=l neZ\{0} ^ ^ 



gfci(si-i)+-+fcd-i(^d-i-i) T{q,s,n,M -1) 

0</si<---<fe<i_2</s<i-i ^ ^ L a ij n£Z\{0} ^ ^ 

/, \ M+1 „oo /I tN — s— M— 1 

•(3) (.s)M+r-/^ ,..».,.(.-2+M) 

^/ci(si-l) + --- + fcd-2(sd-2-l) 

= E rfcj«i...[fc, J«.-2 E W9.fc,5d-i,5). 



where we replace the index kd-i by fc. Similarly, 

Af-l 



C^0=31og,^ Y: , {s)r-SC,{s, 

r=l neZ\{0} ^ ' ^ 



Sd-l + Sd + T"), 



v-^ Cm -!) + ■■ ■ + /.-,, ii.s,, 1-1) r(r7..s-,,)..U- 1) 

i?o--31ogg 2^ [fci^.--[fc,_iP-^ (^^ri^Ip 

0<fci<---<fe<j_2<fcd-i ^ L a ij neZ\{0} 

1 \ M „oc / 1 T \ -s-M 



9-1 



foc /-, _ X \ -s-JVJ 

(s)m • / e2-"-g-(«-2+M) 



and 



r=l neZ\{0} ^ ' ^ ^ 

^2 'sp gfci(si-i)+-+fe<i-i(s<i-i-i) T(g,s,n,M- 1) 

L>i--(logg) 2^ ^ ^ ^ 

0<fci<---<fcd_2<fcd-i ^ ^ L o ij neZ\{0} ^ 

The crucial step next is to control the summations over and show that they converge 

uniformly with respect to q. When < g < 1/2 this is clear. The only non-trivial part is when 
g t 1- So we assume 1/2 < g < 1. Note that 

limT(g, s,n,r) = -— , lim = 1. 

q\\ ' ' ' {2mny 9T1 g - 1 

Lemma 4.1. Let Sd = s = a + iT. Let go = max{l/2, e^'^-^'^)/^} if t > and let go = 1/2 if t < 0. 
Then for all 1 > q > qo and positive integer k we have 



logg 



g-1 



< 2, and \T{q,s,n,r)\ < 



{6ny ■ 



Proof Let /(g) = 2(1 -g) + logg. Then /'(g) = -2 + 1/g < whenever g > go. So /(g) > /(I) = 
whenever 1 > g > go. This implies that 2(1 — g) > — logg whence logg/(g — 1) < 2. 
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To bound T{q, s, n, r) we consider each of its factors in definition H16|l . For each < j < r we 
have 

\2mn+ (s - 1 + j) loggfp = ((cr - l+j)\ogqf + (27rn + rlogg)^ > (27rn + rlogg)^ 

which is independent j. If r < then clearly |27rm + (s — 1 + j) \ogq\ > 6n. If r > then it follows 
from q > e(6-27r)r ^j^^^^ 

2Tm + T log q > 2?™ + 6 — 2n > 6n, 
as desired. □ 

Next we want to bound the integral terms in D-i. Let \ad\ < N and |cr(j-i| < N' for some 
positive integers N and N' . Fix an arbitrary x > k and a positive integer M > 16 + 2N + 
6E-=i'W + l)- Then 



^-fc(M/6-Ar'-l) 



^ \ -S-M-1 



^ qX(M~N-2)-kM/6 



, X M-N+1 

1 - <? 



1 — q J \1 — q^ 

Denote by g{q) the right hand side of the above inequality. 

Lemma 4.2. Let 1/2 < g < 1. Then g{q) is increasing as a function of q so that 

g{q) < \rmg{q) = ^g^^- 

Proof. Taking the logarithmic derivative of g{q) we have 

g'{q) _x{M -N-2)- kM/6 ^ _ ^ _^ ^^xq--\l - g) - (1 - 9^) 



9{q) q (l-g)(l-(?-) 

_{l-q){l- q'=){x{M - N -2)~ kM/6) + {M - N + l){xq='{l - q) - q + g^+i) 
" q{l-q){l-q-) 

whose numerator is denoted by h(q). Then 

h'{q) ^{{x + l)g^ - xg^-i - l)(a;(Af - - 2) - kM/6)) 

+ {M-N + 1) - {x + l)q^) - 1 + {x + 1)9^) . 

Clearly h'{l) = and moreover 

q^-''h"{q) ={x{x + l)q - x{x - 1)) {x{M -N-2)- kM/6)) 

+ {M - N + l)(^x{x{x - 1) - x{x + l)q) + x{x + l)q'j 
^x{x - l)(fcM/6 + 3x) + qx{x + 1){M - N + 1 - {kM/6 + 3x)) 
= {kM/6 + 3x){x^{l -q)- x{l + q)) + qx{x + 1){M - iV + 1) 
>qx{x + 1){M - N + 1) - xil + q){kM/6 + ix) (since l> q) 

q 

>qx^{M/2- N -8) > 



>gx^|Af - iV + 1 - i-t^(3 + M/6)| (since k < x) 



where we used the fact that if q > 1/2 then {1 + q) /q < 3. This implies that h'{q) is increasing so 
that h'{q) < for all 1 > (7 > 1/2 (recall that h'{l) = 0). It follows that h{q) is decreasing. But 
h{l) = so we know h{q) > for all such q. Thus g'{q) > and therefore g{q) is increasing. This 
completes the proof of the lemma. □ 
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We now can bound the innermost sum of Z)_i. From Lemma l2. II Lemma l4.1l and Lemma |4.2I 
we have (if d = 2 then take ko — 1) 



k=l+ka-2 k=l+ka-2 ^ ^ 

, {M + 1)! (M + N\ ^ qk(M/6-N'-l) 



^M-N+l 



E 



M \ M + 1 J ^ k^i^N-N' 

k=l + kd-2 



<(M + l)!(2Af)^^+i J2 1'^ 



_,fe(M/6-JV'-l) 

since 2^'+*^+\(Af + 1) < Att'^Q^'^-^ and M - N > 2. Therefore by Lemma lO 

which converges as proved in Prop. 12.21 

Exactly the same argument appHes to the integral terms in Dq and Di^ which we leave to the 
interested readers. These convergence results imply two things. First we can show by induction 
on d that H15(l gives rise to an analytic continuation of ^^(si, . . . , Sd) as a meromorphic function on 

\ 6'^. Second, also by induction on d, we now can conclude that it's legitimate to take the limit 
q 1 1 inside the sums of Ci, and D^, to get (note that Hmq^i S'^(q{s) = ({s) for any s G C''"^ \ &'d-i 
and any positive integer n) 

1 Is 

\im (q{s ' ,Sd-l,Sd) = -C{s',Sd-l + Sd-l) - -({s',Sd-l + Sd) + — C(s',Sd_i + Sd + l) 

q'\l S — i Z iZ 

M-1 ^ 

T^—yT^i-')r+i-Ci^'^-''d-i + Sd + r + l) 



(27rm 

r=l neZ\{0} ^ 



M+1 



i d-1 „cT\!n\ ^ ' -"Sd-l 



E — f («d)r-l • C(S1, ■ • • , + + r - 1) 



r 

r=0 



1 1 Z""" 5 (■s)m+i 



l,s.,.us,_, / +1 (x) 



(M+1)! ^ K-'-K^^'h. 
by lO and its specialization with x — 1 



E 

nGZ\{0} 



1 _ 5^+2(1) _ Br+2 

(27rm)''+2 ~ ~ (r + 2)! ^ "(r + 2)!' 



The main theorem now mostly follows from Thm. The poles at s^ = to — Tog^"" given by 
the first term in formula (|15|l when m ~ 1 and n — and by the terms T{q, Sd, n, r) as defined in 
(I16|) if m < 1 and n 7^ 0. The location of the other poles are obtained by induction using those 
poles of the g-Riemann zeta function presented in Thm. [TTTl for the initial step. This completes the 
proof of our main theorem. 

5 Series g-shuffle relations 

The classical multiple zeta functions satisfy shuffle relations originating from their series represen- 
tations. For example, 

C(S1)C(S2) = C(S1, S2) + C(S2, Si) + C(S1 + S2). (17) 
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In general we can define a shuffle operation * on finite words so tliat for any complex numbers a, 
b ( regarded as letters) and words wi and W2 of complex numbers 

awi* bw2 = a{wi* bw2) + b[awi* W2) + (a + b)[wi* W2) 

so that treating complex variables as words we have 

C(wi)C(w2) = C{WI*W2). 

To generalize these relations we first define shifting operators on words of complex variables: 
(si, . . . , Sj_i,5(sj), Sj+i, . . . ,Sd) ^ Si{si, . . . ,Sd) 
and set C(5i(si, . . . , Sdf) ~ Si({si, . . . , Sd). Now we define the g-shuffle operator *q on words by 

awi*q bw2 — a{wi*q bw2) + b{awi*q W2) + S{a + b){wi*q W2). 
Theorem 5.1. For any two words of complex variables wi and W2 we have 

Cqiwi)Cq{w2) = Cq{'Wl*q W2) . 

Proof. Induction on the length of wiW2. □ 
For example 

Cg(si)C?(s2) = Cg(si, S2) + Cq(s2, Si) + Cg(si + ^2) + {q - l)Cq{si + ^2 - 1). (18) 
We can recover the shuffle relation of H17|) by taking limq|i. 

6 Special values of Cq{si,S2) 

For integers ni, . . . , rt^ we set 

Cq{ni, . . . ,nd) = lim ■•• lim C« ■ • ■ , Sd), 
C,qini,...,nd)^ lim •■• lim Cq{si, . . . , Sd) 

if the limits exist. Interesting phenomena occur already in the case d — 2 and these should be 
generalized to arbitrary depth. By Thm. we get 

(9(51, S2) -(1-9)'^+"^ 2^ ' 



= (1-9) 



ri J\ r2 y (1 - q'*2+r2-l)(l _ gS2+»-2 + si+ri-2) 

I „2s2+si-3 2S2+S1-2 

S1+S2 I y I £iy 



(1 - (7-''2)(l - q'*2+Sl-l) (1 _ q-->2)(l - gS2+Sl) 

Sl(si + 1)52^2+^1-1 52(S2 + l)(72s2+si + l 



2(1 - g-'2-i)(l - g-'2+si) 2(1 - g''2+i)(i _ 5S2+^i) 
Clearly we have 

1 3 

Cg(0, 0)= lim lim Cq(si,S2) 



' ' .i-Os2^o"«^ (g2_l)(5_l) 2(<7-l)logg log^g' 

Cg'^(0,0)= lim lim Cg(si,S2)=- ^ ^ " 



^2-0.1^0^"^'^^' (g2_i)(g_l) (g-l)logg 2(g-l)logg- 
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It is not too hard to find that 

limCg(0,0) = i hniC,''(0,0) = ^. 

This is consistent with what we found in by using generahzed functions (distributions). Sec 
also equations H22(l below. In ^Jj we further showed that near (0, 0) the double zeta function has 
the following asymptotic expansion: 

/-f \ 4si + 5.S2 

12(si + S2) 

where i?(si,S2) is analytic at (0,0) and lim(5j^ ,^2)^(0,0) R{si,S2) — 0. 

Let n, k be two non-negative integers, and m = k — n — 2. Wc now consider the double zeta 
function around (si, S2) = {^f^, which has the following expression by Thm. 

— ^C(.si + .S2 - 1) - ic(.si + s,) - X: )C(S1 + S2+r-l) 

r=2 ^ ' 

+ a(m)(-l)"|f -n- 2)!(s2 + n)C(s2 + si + A: - 1) (19) 
where aim) = if m < —1 and aim) = 1 if to > 0. Note that the last term is zero when computing 

<^(— TO, — n) = lim lim ('(■511*2) 

si — *— m S2 — >—n 

while it has possibly nontrivial contribution for C^(— to, —n) since 

lim lim (s2 + ")C(*2 + Si + — 1) = 1. 

S2 — ^ — n s\ — ^ — m 

We get 



fc, n 


pole, residue 


indeterminacy, ( = C,^^'{n + 2 — k, —n) 


/c = 


-l/(n+ 1) 


NO 


/c = 1 


-1/2 


NO 


2tfc,3 < fc < n + 1 


NO 




2\k,k = n^2 


NO 


Bfc_i/(fc-l) 


2\k,k> n + 2 


NO 


Bfc_i/2(fc-l) 


2\k,2 <k<n+l 




NO 



Table 1: Poles and indeterminacy of double zeta function. 



When TO > and 2\k the values of C(— to, — n) and C^(— to, — n) are different in general: 

fc(n + 1) r \r — 1/fc — r 

C^(-TO, -n) =C(-"i, -") + i-lT^nUk - n - 2)!. (21) 

Note that the term corresponding to r = 1 is non-zero if and only if fc = 2 (and n = 0). From this 
observation we again recover that 

C(0,0) = i, C^(0,0) = ^. (22) 
We now consider the g-double zeta function. 
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Theorem 6.1. Let k, n be two non-negative integers, and m = k — n ^ 2. If m < —1 then the 
q-double zeta function Cq{si,S2) has a pole at {—m,—n) with residue given by: 



( k 



Res 

(si,S2) = (-m,-n) 



Cq{si,S2) 



(l-g)2-'=(log<z)-i 



= < '-=0 



n + 1 — r \ / n 

k — r 



^, i-=0 



(-1)" 



(n + 1) logg 



if k < n, 
if k = n + I. 



Proof Use Thm. O 

Corollary 6.2. Let n be a non-negative integer. Then 

R-es C(si,S2)^^ -Cq{-n) 

{si,S2) = {l-n) logq 



and 



\imCq{-n) 



Bn+l 



Res C('Si,S2) 



n+1 (si,S2) = (l,-n) ^ 

Proof. Equation H23|) follows from the case to = —1 in the above theorem and |S1 (6)]: 



□ 

(23) 
(24) 



C,(-n) = (1 - q)- J/(9"+'"'" - 1) - 



(-ir 



(n+ l)log9 



The first equality in H24|l is 5, Thm. 1] and the second equality follows from Tabled 
Corollary 6.3. Let k, n be two non-negative integers, m — k — n — 2< —2. Then 
lim Res Cqisi,S2)= Res C(si,S2). 

iTl (si.S2) = {-m,~n) (si,S2) = (-m,-n) 

Proof. By Thm. 16.11 and Table Q we only need to prove 
k 



□ 



(25) 



lim(l-g)i-'=^(-l)' 

r— 



n + 1 — r \ f n 
k — r 



Hi 



,7i-\-l — r 



1) 



-1 



n+l 



if fc = 0, 

if 1 < fc < n. 



(26) 



First by generating function of the Bernoulli numbers 

1 _ _ 

r _ -| p(n+l-r) log g _ -| 

1=0 



/-I 



Plugging this into the left hand side of equation 126|) . replacing 1 — g by — logg, and exchanging 
the summation we get 



(-logg)i-'^^(-l)' 



r=0 



n + 1 — r\ (n 

k — r 



n+l—r 



1) 



(-l)i-'=^f (logg)'-'=^(-l)^ 



r=0 



n + 1 — r \ (n 

k — r 



(n+l-r)'-i. (27) 
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Then the inner sum over r is the coefficient of x*^ of the foUowing polynomial 

n+l k 



k=0 r=0 



J, f n + I — r\ f n 

k-r 



n+l 



n + 1 — r 



(28) 



where y = x/{x + 1). When I = this expression becomes 



[x + 1) 



n+l 



— E(-^r 



n + 1\ {x + 1) 



n+l 



n+l 



n+l 



n+l 



l-(-x) 



n+l 



Note that k < n we see the coefficient of x'^ in /o(a;) is if fc > and it's l/{n + 1) if fc 0. If fc = 
then only the constant term — l/(n + 1) in H27|l remains when q ] 1 which proves the corollary in 
this case. So we can assume l,k > 0. Then 



= (.^) " {(. + 1)-^ ±i-yr (;:)-"+^-^} 

-1 

[{x + iy^+'ziz-yr} 



z = \ 



Note that highest degree term in ]i{x) is contained in 



\n-Z+l 



=n(n - 1) • ■ • (n - ? + 2)(x + 1)"+'(1 - y) 
^n(n - 1) • • • (n - ? + 2)(.t + 1)'. 

If ^ = 1 one can easily modify this to get just x + 1. If / < fc then the coefficient of x^ in ]i{x) is 0. 
If ^ = fc it is equal to 

n 

-fc- ly 



i(n - 1) • • • (n - fc + 2) = (fc - 1)! 



The last express is valid even for fc = ^ = 1. Thus the range of I in the outer sum of (|27|1 starts 
from fc. Moreover, the first term of (|27|l is 



as desired. This completes the proof of the corollary. 



□ 



Proposition 6.4. Let fc, n he two non-negative integers such that n > k and fc is even. Let 
m ~ k — n — 2. Then 



Res Cq{si,S2) 
(sx,S2) = (-m,-n) 



-/(g)(g-l)/logg 
D{q) 



n+l 



J— n+l — /c 



where f{q) G 1j[q] is a palindrome with leading coefficient (^), F{q,j) G Z[q] is a factor of (q-' 
l)/{q — 1), e — or 1, and deg^ D{q) = n + deg^ f{q), such that 

lim Res Cqisi,S2)= Res C(si)S2)- 

qll {si,S2) = {~rn, — n) (si ,S2 ) = ( — m,— n) 
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Proof. The computational proof is left as an exercise for the interested readers. □ 
Example 6.5. By Thm. IHTI wc find with the help of Maple 

(3„.S,-4^^'^''^^ - Pi(<z,2)Pi(g,3)Fi(g,4) " 
where Pa{q, m) = X^JLo Moreover we can check that 

lim Res C<3(si,S2)= Res C(si: sa) = --^ 

by Table n with k = 2 and n = 4. 
Example 6.6. From Thm. IHTI we get 

-Uq^g{q){q- I)/ log q 



Res Cqisi,S2) 



isus2)He.-sr" ' Pi(g,4)Pi(q,5)Fi(g,6)P2(g,3)P3(g,2) 
where gijq) is a polynomial in q of degree 14 satisfying 

/(g) = 5(7" + 6(7!^ + 8q^2 _^ 7^11 _ ^10 _ 20g9 _ 30g8 - 34g^ 
Then we can compute with Maple 

7 

lim Res ((si,S2)= Res C(si:S2) = -r 

■ZTI (si,S2) = (6,-8) ^ (si,S2) = (6,-8) 15 

by TableHwith A: = 4 and n 8. 

Example 6.7. Consider the point (si,S2) = (5, —9). We have 

-\2qS{q){q-\)l\ogq 



Res C,q{s\,S2) 



(.i,.2T=('5,-9)'''"""'' Pi(g,4)Pi(g,6)Pi(g,7)P2(<Z,2)F3(<Z,2)A(g,4) 
where m) = X^j'Lol^l )"''?"' ^-iid g{ci) is a polynomial in q of degree 18 satisfying 

q^^g{\lq) = giq) ^ 2q^^ - - Iq^'" - llg" - 16g" - Aq^^ + 9g" + 28^^" + 30g^ + • • • 
so that we again have the equality 

lim Res Cq{si,S2)= Res C(si- ^a) = 

9TI (si.S2) = (5 -9) ^ (si,S2) = (5 -9) 2 

by Table^with k = 6 and n = 9. 

Theorem 6.8. Let m, n &e two non-negative integers and k — m + n + 2. Then Cglsii^a) has 
indeterminacy at (— m, — n) swc/i t/iat 



Cq{-m,-n) =(1 -g) 



(-1)'= 



(to+ l)(n + l)(logg)2 



r=0 



{n + 1) log g V ^ / g'"+i"'' - 1 



« (-1)'-+™+! m!(n + l-r)\ fn\ 1 



T'=0 

m n 



logg (fc — r)! \r J q"^+^~^ — 1 

1 1 



n— r2— 



14 



and 



+ X] /^fc - n - 2^^ n!(m + 1 - r)! 



r=0 
m n 



logq \ r J (fc-r)! q-m+i-r _ i 

1 1 



fi/ \r2j - 1 qk-ri-r2 _ i 



n— r2— 

Proof. Use Thm. ITl 

Similar to Cor. 16.21 and Cor. 16. 31 we have 
Corollary 6.9. Let m and n be two non-negative integers. Then 

limCq(-m, -n) = ((51,^2), limCg^(-m, -n) = C^(si,S2)- 

Proof. Set k = m + n + 2. We consider Qq^{—m, —n) first. From Thm. I^TSl we get 



* vogq 



9-1 



2 



where 



00 m , ^-.^ 



i=0 r=0 

BjB 



We first compute B as follows. Write 



4=0 

where 

If i = then we can prove by decreasing induction on n that 

i!n! 1 



W(m,n,0) = ^(-l)'■+''■+"- 



r!(fc-r)! fc(n+l)' 

This is trivial if n = fc — 2. Suppose H30|l is true for n > 1 then we have 

k — n — 1 (n — lY 1 

W{m, n - 1, 0) = W{m, n, 0) + , ^ ! = — , 

n [n + 1)\ kn 

as desired. 

Similarly we can compute C as follows. Put 

c - 1 t ± (:;) c: ^^)-(^ - 
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Wc now change the upper hmit of r2 from n to n + 1 in the above. The extra terms correspond to 
those by setting i = and r2 = n + 1, which produce exactly A. Therefore, 



B^B, (log q) 



i+j 



-V{m,n,i,j) ~ A 



where 



n— r2— 



n + 1 \ri 



r2 



(n + 1 - r2)''(fc - ri - ra) 



For J > 1 we have 



V{m,n,iJ) = X 



dx 



dyj 



(-1) 



fe+ri+r2 



n— r2— 

to\ /n + I 
ri) \ r2 



n + 1 



„,n+l-i"2 fc-r-i -r2 



(31) 





(-l)'=n! 
(n + 1 - «)! 



it i + j < k, 
{k — i — 1)\ it i + j = k,i < n + 1. 



We have used the fact that if i + j = k and i > n + 1 then j < k — n — 1 and by exchanging 
the two operators x{d/dx) and y{d/dy) we can easily show that (|31f) is zero. So if Z < fc the total 
contribution to the coefhcient of (logg)' from l/(m, n, j) with j > is trivial and if / = fc it is 
equal to 



Bk-i 



n+l 



BiBk 



i=0 



(-l)''n!(fc-^-l)! _ J 2(/c - 1) 
i\{k - i)\{n + 1 - i)l " 



Bh 



ri+l 

E 



BiBk-i I n 



k{n + 1) -^-^ i(fc — i) \i — 1 



if 2 t fc, 



if 2 I fc, 



(32) 



because k > n + 2 > 2 and Bk — it k is odd. 

To deal with V{m, n, i, 0) note that equation H31|) still makes sense if we interpret the operator 
x{d/dx)^^ as follows: 



4) {FM}\ ^fm^ 

ax J ^ J x=i Jq X 



whenever F(0) — 0. Thus we get 



V{m,n,i,0)^^-^(y4- 



n+l \ dy 



(xy- l)"+^(x- l)™da; 



Therefore if i = then we get 



y(m,n, 0,0) = 



fc rl 



(-1) 

n+l Jq 



{x~l)''-Ux 



1 



A;(n + 1) 



= -W{m,n,0) 



(33) 
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from equation (|30|l . If i > 1 then integrating by parts we get 



y{n + 2) 



y{n + 2) Jo 



{xy~lY+\x-ir-'dx 



= (-1) 



= (-1) 



fc+i 



y{n + 2) y'^{n + 2){n + 'i) 



+ --- + (-ir 



TO!(n + 1)! 
y™(TO + 71 + 1)! 



m\{n + 1)! 



+ (-1)' 
^ ' y"^{m + n+l)\ jQ 

fc+„m!(n + l)!(y-l)fc 



(x-y - 1) 



,m+l 



fc! -y 

It follows from changing the index r to m — r that 



m!(n + 1)! 



y''+^{m - r)!(n + 2 + r)! ' 



F(m,n,i,0) ==(-1) 



fe+ri ' 



A:! \^ dy 



y 



m+l 



2/=l r=0 



r\{k — r)\ 



(m + l - rf 



{-l)''+^W{m,n,0) ifO <i<k, 

{-l)''+^W{m, n, 0) + (-l)'^+"n!(fc - n - 2)! if i = A:. 



(34) 



Thus when < i < k and k is even we have V{m, n, i, 0) = —W{m, n, i). It follows from H32|) . (|33|l 
and (E3)) that 

limCg^(-"l, -n) = C^('Sl7S2) 

since = if fc > 2 is odd. 

Let's turn to prove the first equahty in if^ . By Thm. l^^Hl we have 



Cq{-m,-n) ^{q- ly 



\ogq 



{D + A + E + C) 



where A and C are as above and 



D 



where 



Hence 



(m + l)(n+ 1) 

n 

C/(m,n,z) = -^(-l) 
[/(m,n,0) = T:^(n,m,0) 



m, n, ij, 



i=0 



r+n 



7,lr7,' 



r=0 



"[(fc -r)! 



(n + 1 -r)\ 



-1 



1 



fc(m+l) fc(fc-n-l) fc(n + l) 



D. 



We only need to show that 



U{m, 71, i) 



if < 7 < fc, 



W{m, 71, z) 

W^(77i,7i,7) - (-l)"n!(fc - 71 - 2)! ifi = fc. 



(35) 



(36) 
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Indeed when i > we have 



,n+l — r 



dv J I fc! V?" 



n+l — r 



m!n! , ^ 

-{n+l-ry. 



r!(fc — r)! 



Then first term is if i < A; and it's (— l)"+^m!n! if i = fc. When r = n + 1 the summand in the 
second term is zero since z > 0. So we can let r range only from n + 2 to fc. Then change the index 
r to fc — r (and let r run from to m) we can see immediately that the second term is exactly 
W{m,n,i). This proves equation (|36() which together with (|35() implies the first equation in H29(l . 
We thus finish the proof of the corollary. □ 

We conclude this section by remarking that by shuffle relation (|18|l we can also analyze Cq{si,S2) 
at {—n, n + 2 — k) for any non-negative integers fc and n. For example, it's easy to compute directly 
that 

1 q-^+l 



Res Q(si,S2)=-j. \-, E(-1)'PV^ + 1)t 

,2)={-3.2)^' (l-g)logg^' Vv 1 

-qiq-l 



{q + l){q^ + l){q^ + q + l)\ogq' 
which can be obtained also by the shuffle relation (|18|) and the expression 

qiq-lf 



Res Cq{si,S2) 



(si,..)=(2,-3r^ ^' iq+l){q^ + l)iq^ + q + l)logq 

by taking fc = ri = 3 in Thm. Ki.ll Thus (—3, 2) is a simple pole of Cg(si: sa)- On the other hand 
C(si, S2) does not have a pole along si + S2 = —1. Indeed we find that 

lim Res {C<j(si: ^2)] = 0. 

7 g-multiple polylogarithms 

It is well known that special values of the multiple zeta function ('(■Si: ■ • ■ i Sd) at positive integers 
(ni, . . . , nd) can be regarded as single- valued version of multiple polylogarithm ii„j...._„^(zi, . . . , Zd) 
evaluated a,t zi = ■ ■ ■ = Zd = 1. For \zj \ < 1 these functions can be defined as 



Lini,...,na{^lT ■ ■ ■ 1 ^d) — ^ ^ 

By Chen's iterated integral 



•^1 • ■ ■ 



fc"^ fc"'' ' 
0<fei<-<d<j 1 ■■■ d 



Jo h - ai \ h J td- ad \td J 

where aj — \/ JliLj -^j f*^^ j = 1 , . • . , d, we can obtain the analytic continuation of this function 
as a multi-valued function on \ where 
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When \zj\ < we define its g-analog (0 < g < 1) by 



•^1 ■ ■ ■ 



0<fci<---<di ^ ^ ''J 

Clearly when (7 | 1 we recover the ordinary multiple polylogarithm. Moreover, the special value of 
g-multiple zeta function at positive integers 

Cq(ni, . . . , nd) = iig;„i, • • • , g"""^)- 

Note that our definition of the g-multiple polylogarithms is different from that of 8 . In case of 
logarithm and dilogarithm our definitions are different from that of T. We want to convince the 
readers that ours are also good analogs of the ordinary ones. 

We can mimic the method in section|5]to get the analytic continuation of iig;ni,....nd (^ii ■ ■ ■ 1 ^d)- 

Theorem 7.1. The q-multiple polylogarithm function Liq-„-^^.. .,na{zi, ■ ■ ■ , Zd) converges if \zj\ < 1 
for all j ^ 1, . . . ,d. It can be analytically continued to a multi-valued function over \ Sg^^ via 
the series expansion 



+OC d 

Lz,^„,,...,„,(zi,...,z,) = (l-g)«^+-+"^ ^ J] 

ri,...,rd=0 j=l 



(38) 



Proof. The first part of the lemma is obvious. Let's concentrate on the analytic continuation. By 
binomial expansion (1 — = X]^o {^'^^r~^)^^ 



+00 d / 

+ r, — 1 



i%.......„.(zi,...,..) = (i-gr+-+- E E n ' Y^^'f'f' 

0<fei<---<fc<j ri,...,rd=Oi=l ^ ^ ^ 

As < (7 < 1 the series converges absolutely by Stirling's formula so we can exchange the summa- 
tions. The theorem follows immediately from Lemma 12.31 bv taking Xj — Zjq^^ . □ 

However, this analytic continuation is not suitable for comparing with its ordinary counterpart. 
The we define Jackson's g-diffcrential operator (cf. (31) by 

(1 - q)z 

Lemma 7.2. Let d,ni, . . . ,nd be positive integers. If nj > 2 then we have 

I^q;Zj Liqyrii ,...,71(^(^1; . . . i — Liq^rii ,. . .,nj — 1,. . .jn^j (^^1 ; . . . 7 Z(l) j 

^3 

if d > 2 and nj — 1 then 

Dq;ZjL'iq;ni,...,na{zi, . . . , Zd) — - -^^g;ni,...,ri] ,...,nd (^1 , . . . , Zj-lZj, . . . , Zd) 

i — Zj 

T^j r-^^g;ni,...,nj ,...,nd (-^1 ; ■ ■ ■ ^ ZjZj^i^ . . . ^ Zd\ 

Here the second term does not appear if j = d. If d — ni = \ then 

Dq.^Liq.liz) = 

Proof Clear. □ 
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The same properties listed in the lemma are satisfied by the ordinary multiple poly logarithms. 
We note that the first equation in 8, Lemma 1] is valid only for nj > 2. 

Recall that for any continuous function f{x) on [a,b] Jackson's g-integral (cf. ^) is defined by 

/ fix) d,x := /(« + a)){q' - q'+'Kb - a). 

Then for every x G [a,b] we have 

/ D,,tf{t)dgt = fix)- f{a), and lim / f{t)d,t=[ f{i)dt. 

Remark 7.3. Note that in general f{x) dqX + f{x) dqX ^ f{x) dqX. 

Similar to Chen's iterated integrals we can define the g-iterated integrals as follows: 



(39) 



J a ti - ai 



dqU _ ( ( ['^ dqh \ dqt2 dgU-l \ dgU 



tr — ar J a \J a J a \J a ti — Gi J t2 — 02 tr-1 — O.r-1 J tr — Ur 

Define 

d 

^q;d =■■ {{zi, ...,Zd) eU' ■.Y[z, = q-"^, m G Z>o, J = 1, . . . , d}. 

Corollary 7.4. We can analytically continue Liq.ni,...^nd{zi, ■ ■ ■ , ^d) io \ 'Dq-^ by the q-iterated 
integral 



i- 



-^^d d,qtl ^ / dqtl 

Jo ti - ai \ ti 



dqt, fdqt,^°^^'^-'^ 
o . . . o o ' 



td — ad V id 



where aj ~ 1/ Jlf=j — Ij ■ • ■ ? d,. Further, for all (zi, . . . , Zd) G such that Jl^Lj Zi ^ 

[1, +oo) for 1 < j < d we have 

lim Liq-n^ rid {zi, ■ ■ ■ , Zd) = Li^ rid {zi, ■ ■ ■ , Zd), 

where Lin-^....^ndizi, . . . ,Zd) is defined by ()37|l with the path being the straight line segment from 
to 1 in . 

Proof. It follows from Lemma 17.21 and equation (|39|l . The singular set 'Z)q-^d is determined by 
Thm. mi so that for each j = 1, . . . ,d the function l/{t — aj) is continuous on [0, 1]. □ 



8 Iterated integral g-shuffle relations 

In section|Slwe encountered some g-shuffle relations of g-multiple zeta functions. Classically, multiple 
zeta values satisfy another kind of shuffle relation coming from their representations by Chen's 
iterated integrals. In our setting we have seen that special values of q-multiple zeta functions can 
be also represented by g-iterated integrals. In this last section we would like to study the shuffle 
relations related to these g-iterated integrals. We shall see that they're more involved than their 
ordinary counterparts. We start by writing 

Shfl(ui 0---OUr, Ur+1 O • • • O Ur+s) = ^ W<t(1) O • • • O U„^r+s)- 

a 

where a runs through all the permutations of {1, . . . , r + s} such that a~^{a) < a~^{b) whenever 
l<a<b<roTr + l<a<b<r + s. For any expressions Fi we put 

r 

y = Fi o . . . o 
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Lemma 8.1. Let Ui = dqt/{t — ai) and Vj = dqt/[t — bt) where \ai\, \bj\ < 1 for all 1 < i < r and 
1 < J < s- Let a be any positive number. Then 

na j-a j-a min(r,s) 

/ Ui o ■ • ■ o ■ / vio---oVs— I Sh&[ui o ■ ■ ■ o Ur,Vi o ■ ■ ■ o Vg) + [q — iy- 

Jo Jo Jo 

o • • • o Ui^-i,vi+j^_^ o • • • o Vj^-i) o {ui^,Vj^), (40) 

<-<ic<'- ° a=l 



l<il<-<ic 

l<jl<---<Ja<S 



where io = jo = 0, v+i =r + 1, jc+i = s + 1, {ur+i,Vs+i) = 1, and for all i,j 



Proof. This can be proved by induction on r + s. We only want to mention that the key formula is 
Dg.Afix)g{x)] = [Dq.,fixMx)+fix)[Dq,,gix)]+x{q - l)[Dq,,fix)][Dq,,gix)]. 

□ 

We will say the term {ui, Vj) is a collapse in the shuffle. The lemma roughly says that iterated 
shuffle relations is different from those produced by Chen's iterated integrals because collapses may 
occur. The number of collapses is at most min(r, s). 

Lemma l8 . 1 1 implies that if m, n > 2 are different then 



1 d„t rdgty^""-'^ dgt (dgtv^'^-^'^ 



Ag{ni,n) + Aq{n,m) + Bq{m,n), (41) 



where 



m-lmin(a,n) i / j , \ o(m-l-a) , , /j . \ o(n+a-l-c) 

v-^ v-^ , , / dnt / d„t\ ^ ' d„t / d„t\ ^ ' 



± .... .. ^ 

Aq{m,n) = Y ^ E{a,n;c) 

a=0 c=0 



t-q^-"^ V ^ / t-q^-"^ V ^ 



Bq{m,n) —{q — 1) E{m — l,n — l;c) / 

c=o -^0 



(t-gi-")(i V * 

Here the coefficient E{r, s; c) represents {q — ly times the numbers of terms in shuffle of Ui o • ■ • o 
and vi o ■ ■ ■ o Vs with c collapses (see (|40ll ) . It is not hard to see that E{r, s; 0) = {^'l'^)- Thus 



E{r,s;c)^iq-ir ^ 11 



C+l . _ . _ . _ 



l<il< - <ic<i' a=l 
l<jl<-<Jc<s 



ia ia—1 1 



We want to convert the expressions in H41() into something that is close to linear combinations of 
multiple zeta functions. Since m ^ n we get 

tdqt 



(i-9i-™)(t-9i-") V t 

1 dqt 1 dqt \ 



k=l ^ ' ^ k=l ^ J 
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Proposition 8.2. For any positive integers m, n we have 



min(m,n) — 1 

Bq{m,n) ~ {q ~ I) ^ E{m ~ l,n - l;c)- 



c=0 

1 

r'?7i — 1 — c 



( -5"-i-^Cg(TO + n - 1 - c) + 



Proof. It follows from Cor. 

To handle Aq(rn,n) we need to evaluate 



Liq-m-a,l3{q"^ 



,^ (it ^+ /rfA°(^~^^ 



O — — O — O 



/o t-q^-""- V ^ / t-q^-"^ \ t 

where < a < m — 1 and max(n, a) < /3 < n + a. By Cor. 12. 61 we get 

iz,;™-a./3(g™-", g"-^) = T (1 - qyLtq.,,n-c.,p^:{q"''\ q^^'~'). 

So we need to evaluate 

„im-n)k J-f-l)l 

Li fo™"" cP^^) - \ ^ ^ 

l^lq-ra-a,-iyq ,q ) - [/t]™-" ' 



If n > a then by Cor. 12.61 we get 

Liq,^-aMr-'',q''-')= E r " )il-qyC,{m-a-j,j). 

j=0 V / 

By taking a = a, /3 = n + a — c and 7 = /? — i in 14211 and (|43|l we get 
Proposition 8.3. If n > rn then 

a=0 c=0 1=0 j=0 V * / \ 3 J 

We now consider the case n < m. For j > define 
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Set rOCq(7) = limj^o^^Cg(7) = E/>i(^ " l)q'^'^'^'>' /[l]'^ . To evaluate these we need 
Lemma 8.4. For any e > we have 

1>1 ^ ' ! = ^ ^ i>l 1^ J 

Proof. Set a(e,7) = Liq.^{q'^^'^). If e = we have 

a(0, 7) = = - 1)^C,(7 - ^) + (g - 1)''-' E I' 

1=0 01 ^ ^ 
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which confirms the lemma in this case. If e > 1 we have 

T-S Je+l)l 

a(e, 7) = E (9 - -hl~j) + iq- 1)^"^ E 

i=o i>i ^ ^ 



j=0 4=0 ^ ^ i>l 4=0 



7-2 



(T=0 ^ ^ />1 ^ J 



by induction and the combinatorial identity 



4=0 



Corollary 8.5. For all j > 1 we have 

Now we fix 7 and for all s > r > 1 we set 

X^{r, s) := X{r, s) Liq-mil'' , q'^'^)- 
Lemma 8.6. For every e > we have 

X,{r, r + e) = Y^{q If (' + C,{r -^,^) + {q^ If E 

4=0 ^ ^ j = 

Proof. Clearly 



X(r,r) =C?(^,7)+ E 



l<k<l ^ ' ^ ' 

=Ur,l) + iq-l)X{r-l,r-l) 

- E('^ - - ^, 7) + - l)-^ E mf"' - 

4=0 ;>i ^ ^ 

And for e > 1 



X{r, r + e) ^X{r, r + e - 1) + (q - l)X{r - 1, r + e - 1) 

r-l 

=E(9-ir^(^-^^+e-i-*) + ('z-ir E 



(7-1); 



4=0 i<fc<; ^ J 

'"-1 l^e _ ^el qi"f-l)l 



□ 



= E(9 - -i,r + e-l-i) + {q- ly E Y 

4=0 1>1 

=T.(i-^y\ E (^^!7^)(9-i)%(^-^-j,7)| +(9-irE^'c,(7) 

4=0 y j=o ^ ^ ^ J i=o 

by induction. Let a — i + j in the first sum of the last expression. Then the lemma follows from 
(113. □ 
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Define = J2i>i q^^+^'>^ /[I? ■ We now have 
Proposition 8.7. If n < m then 

n — 1 min(a,n) a — c n—a — 1 



Aq{m,n)^Yl E E E ^(a:";c)r . j Cg(m - a - j, n + a - c - i) 

^— n — n 1 — n — n V ^ / \ 3 / 



a=0 c=0 1=0 j=0 
m-1 min(a,ri) a-c 

+ ^£'(a,n;c)( . \{l-qyXn+a-c-i{m-a,m-n), 



a—n c— i— 



+ (9 - !)'■ E ^^^J^fe(j) - C.(2)) - (9 - 1)^(5 - r + 1)^(7). 

Proof. This is clear from equations (|42|l and H43(l . Lemma [8.61 and Cor. 18.51 □ 

Putting everything together we arrive at 

Theorem 8.8. Let m ^ n be two positive integers no less than 2. Then 

(q{m)Cq{n) = Aq{m,n) + Aq{n,m) + Bq{m,n), 

where Aq{m, n) is given by Prop. and Prop. |<y. 7| and Bq{ni, n) is given by Prop. \8.HA 

It is not hard to see that when g t 1 we recover the ordinary shuffle relations of the MZVs 
originally produced by using Chen's iterated integrals. The only unpleasant terms in our q-analogs 
are given by £,q{j) which is closely related to C,q{2). 
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